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Abstract 

The folding of the triangular lattice embedded in two dimensions (discrete planar folding) is 
investigated numerically. As the bending rigidity K varies, the planar folding exhibits a series 

<p : 

of crumpling transitions at K ~ —0.3 and K ~ 0.1. By means of the transfer-matrix method 

for the system sizes L < 14, we analyze the singularity of the transition at K ~ —0.3. As a 

result, we estimate the transition point and the latent heat as K = —0.270(2) and Q = 0.043(10), 
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respectively This result suggests that the singularity belongs to a weak-first-order transition. 



At sufficiently low temperatures, a polymerized membrane becomes flattened macroscop- 
ically [lj; see Refs. [2|-l5| for a review. It still remains unclear Ja-81 whether the crumpling 
transition (separating the flat and crumpled phases) is critical |9|-|21| or belongs to a discon- 
tinuous one with an appreciable latent heat 22N24|. 

In this paper, we investigate a discretized version of the polymerized membrane embedded 
in two dimensions 25H28|; details are overviewed afterward. This model, the so-called 
discrete planar folding, exhibits a series of crumpling transitions at K m —0.3 and 0.1 



123, 



as the bending rigidity K changes. The latter transition exhibits a pronounced 
discontinuous character, whereas the nature of the former transition remains unclear. In 
this paper, we utilized the transfer-matrix method 27| for the system sizes L < 14. We 
implemented a modified folding rule 29], Eq. (J5J), which enables us to impose the periodic- 
boundary condition. Technically, the restoration of the translational symmetry admits a 
substantial reduction of the transfer-matrix size. 



To begin with, we explain a basic feature of the discrete planar folding [27|, |28|; see 
Fig. [1] (a). We consider a sheet of the triangular lattice. Along the edges, the sheet 
folds up. The fold angle 9 is either 6 = (complete fold) or it (no fold). The elastic 
energy at each edge is given by KcosO with the bending rigidity K. The thermodynamic 
property of the planar folding has been studied extensively 27|, |28|]. The transfer-matrix 
simulation for the system sizes L < 9 [27J revealed a series of crumpling transitions at 
K rj —0.3 and K = 0.11(1). The behavior of the specific heat around K ~ —0.3 indicates 
that this transition would be a continuous one. The cluster variation method (CVM) of a 
single-hexagon-cluster approximation 28| indicates that there occur crumpling transitions 
at K = —0.284 and K = 0.1013 of the continuous and discontinuous characters, respectively. 

The crumpling transition K w —0.3 is closely related 32J to that of an extended folding 
30h32| at K 3 sa —0.8. (The extended folding, the so-called three-dimensional folding, has 



four possibilities, cos# = ±1, ±1/3, as to the joint angle 9.) That is, according to an 
argument based on a truncation of the configuration space 32|, the following (approximate) 
relations should hold; 



K = K 3 /3 

Q = Q 3 - 



(1) 

(2) 



Here, the variables Q and Q3 denote the latent heat for the planar- and three-dimensional- 



folding models, respectively. A number of results, (K 3 ,Q 3 ) = (-0.852,0) [32], (-0.76(1), 0.03(2)) 



33| . and (—0.76(10), 0.05(5)) J29j, have been obtained via the CVM, density-matrix renor- 
malization group, and exact-diagonalization analyses, respectively. The nature of its tran- 
sition at K 3 ps —0.8 is not fully clarified, because the three-dimensional folding is compu- 
tationally demanding. It is a purpose of this paper to shed light on this longstanding issue 
from the viewpoint of the planar folding. (It has to be mentioned that the planar folding 



has to a relevance to a wide class of systems 26|, |34j-|37j . ) 



For the sake of selfconsistency, we present the transfer-matrix formalism for the discrete 
planar folding explicitly. We place the Ising variables {a{\ at each triangle % (rather than each 
joint); see Fig. [I] (a). Hereafter, we consider the spin model on the dual (hexagonal) lattice. 
The Ising-spin configuration specifies each joint angle between the adjacent triangles. That 
is, provided that the spins are (anti)parallel, <Ji<jj = 1 (—1), for a pair of adjacent neighbors, 
i and j, the joint angle is 9 = tt (0). The spin configuration is subjected to a constraint 
(folding rule); the pref actor of the transfer-matrix element, Eq. fl3]), enforces the constraint. 
As a consequence, the discrete folding reduces to an Ising model on the hexagonal lattice. 
In Fig. [U (b), a drawing of the transfer-matrix strip is presented. The row-to- row statistical 
weight Tjo-jjo-'} yields the transfer-matrix element. The transfer-matrix element for the strip 
length L is given by [27J 

L L 

T W t },{<Ti} = (n.8(v2i + V2i+l+(T2i+2 + (T2i-l + V2i + V2i+l mod 3 > °)) ex P(~ J2 H i( K )/ T )> ( 3 ) 
i=\ i=l 

with the local Hamiltonian 
K 

H i(fy = -(<?2iCT2i+l + 0"2i+l°2i+2 + 0"2i+20" 2 i+l + a 2i+l a 2i + <^2i (X 2i-\ + ^i-l^i), ( 4 ) 

due to the b ending-energy cost for spins surrounding each hexagon i. Here, the parameter 
K denotes the bending rigidity, and the expression 5(n, m) is Kronecker's symbol. The 
periodic-boundary condition cr^+j = <Tj is imposed. We set T = 1, considering it as a unit of 
energy. 

In practice, the above scheme does not work. The folding rule is too restrictive to im- 
pose the periodic-boundary condition. So far, the open-boundary condition has been imple- 
mented; more specifically, the ran ge o f the running index % in Eq. ([3]) was set to 1 < i < L — 1 



27| . In this paper, following Ref. 29|, we make a modification as to the constraint (prefactor 



of Eq. (J3J)) to surmount the difficulty. We replace the above expression with 

1 L 
T KI-K} = T YKU. S{a2i+a2i + i+a 2i +2+(T / 2i-i+(T2i+(T / 2i + i mod 3, 0)) exp(- £ Hi(K)-Hi(K')). 

*-* 1=1 i+l ijtl 

(5) 
That is, the constraint is released at a defect hexagon i — I. Additionally, the local bending 
rigidity at the defect is set to K'. In order to improve the finite-size behavior, we adjust K' 
to 

K' = 2K. (6) 

A justification is shown afterward. 

Based on the transfer- matrix formalism with a modified folding rule (jSJ), we simulated the 
planar folding numerically. The numerical diagonalization was performed within a subspace 
specified by the wave number k = and the parity even; here, we made use of the spin- 
inversion symmetry o~j — > — 0$. 

In Fig. [21 we plot the free-energy gap 

A/ = f 2 - fu (7) 

for the bending rigidity K and various system sizes L = 6, 7, . . . , 14. Here, the free energy 
per unit cell is given by fc = —haA i /(2L) with the (sub)dominant eigenvalue Ai( 2 ) of the 
transfer matrix. [Here, the unit cell stands for a triangle of the original lattice rather than 
a hexagon of the dual lattice; see Fig. [TJ] From Fig. El we see a signature of a crumpling 



transition (closure of A/) at K w —0.27. The 



consistent with the preceding estimates [271. 1281] . 



ocation of the transition point appears to be 



In Fig. [31 the approximate transition point K(L) is plotted for 1/L 2 and 6 < L < 14. 
The approximate transition point minimizes A/; namely the relation 

d K Af\ K=K{L) = 0, (8) 

holds. The least-squares fit to a series of results for L = 6, 9, 12 yields an estimate K = 
—0.2697(12) in the thermodynamic limit L — > oo. Similarly, for L — 1,2 mod 3, we obtain 
K = —0.2695(14). (An observation that the data L = and 1,2 mod 3 behave differently 



was noted in Ref. [27[.) The above independent results appear to be consistent with each 
other, validating the l/L 2 -extrapolation scheme. As a result, we estimate the transition 
point as 

K = -0.270(2). (9) 



We then proceed to estimate the amount of the latent heat with Hamer's method [38] . A 
basis of this method is as follows. At the first-order transition point, the low- lying spectrum 
of the transfer matrix exhibits a level crossing, and the discontinuity (sudden drop) of the 
slope reflects a release of the latent heat. However, the finite-size artifact (level repulsion) 



smears out the singularity. According to Hamer [38[ , regarding the low- lying levels as nearly 
degenerate, one can resort to the perturbation theory of the degenerated case, and calculate 
the level splitting (discontinuity of slope) explicitly. To be specific, we consider the matrix 



v= ' , (10) 




with Vij = (i\dj(T\j) and the transfer matrix T. The bases |1) and |2) are the (nearly 
degenerate) eigenvectors of T with the eigenvalues A 12 , respectively. The states {\i)} are 
normalized so as to satisfy (i\T\i) = 1. According to the perturbation theory, the eigenvalues 
of Eq. (JTU1) yield the level-splitting slopes due to K. Hence, the latent heat (per unit cell) 
is given by a product of this discontinuity and the coupling constant K(L) 

Q(L) = \K{L)\y/(V n ~ V 22 ) 2 + W l2 V 2l ^ (11) 

for the system size L. 

In Fig. HI we plot the latent heat Q ( TTTjl for 1/L 2 and 6 < L < 14. The least-squares 
fit for L = 6,9, 12 yields an estimate Q = 0.0482(59) in the thermodynamic limit L —¥ oo. 
Similarly, for L = 1, 2 mod 3, we obtain Q = 0.0391(38). Considering the deviation of these 
results as a possible systematic error, we obtain 

Q = 0.043(10). (12) 

The error margin covers both the statistical and systematic errors. 

We consider the 1 //^-extrapolation scheme. The finite-size data are expected to converge 
rapidly (exponentially) to the thermodynamic limit around the first-order transition point 
for periodic boundary conditions, because the correlation length (typical length scale) £ 
remains finite. Hence, the dominant finite-size corrections in our case should be described 
by 1/L 2 (rather than 1/L). On the one hand, the curve of Fig. HJappears to be concave down, 
indicating an existence of a correction of 0(1 /L). However, this possibility (second-order 
phase transition) should be excluded: In a preliminary stage, we made a finite-size-scaling 



analysis, and arrived at a conclusion that the scaling theory does not apply; the critical 
index v estimated from the excitation gap tends to diverge as L — > oo. Therefore, we set 
the abscissa scale of Fig. @] to 1/L 2 ; actually, the result of Fig. [3] demonstrates that the 
abscissa scale 1/L 2 is sensible. 

As a comparison, we provide a simulation result, setting the defect parameter to K' = 
tentatively. In Fig. [5j we present the free-energy gap Af for the bending rigidity K; the 
scale of K is the same as that of Fig. [2j Clearly, the data of Fig. [5] are less conclusive. As 
a matter of fact, the signatures of the crumpling transition strongly depend on the system 
size L. This result indicates that the choice of the defect parameter K' affects the finite-size 
behavior. In the preliminary stage, we survey a parameter space of K', and arrive at a 
conclusion that the above choice, Eq. (jBJ), is an optimal one. 

In summary, the crumpling transition of the discrete planar folding in the K < regime 
was investigated with the transfer-matrix method for L < 14. We adopted a modified-folding 
rule (jSJ), which enables us to implement the periodic-boundary condition. As a result, we 
estimate the transition point and the latent heat as K = —0.270(2) and Q = 0.043(10), 
respectively. The planar- and three-dimensional-folding models are closely related; see Eqs. 



flTJ and (J2J). Making use of K 3 = —0.76(1) 33] and the present result K = 0.270(2), we arrive 



at K 3 /K = 2.815(43) (~ 3). The relation (pQ) appears to hold satisfactorily; a slight deviation 
indicates that the truncation of the configuration space is not exactly validated. Encouraged 
by this result, we estimate Q3 = 0.043(10) via Eq. (j2j). This result is consistent with 
Q 3 = 0.03(2) [33J and Q 3 = 0.05(5) [29], indicating that the singularity belongs to a weak- 
first-order transition rather definitely. Because a direct approach to the three-dimensional 
folding is computationally demanding, an indirect information from the planar folding would 
be valuable. A further justification of the configuration-space truncation would be desirable 
to confirm this claim. This problem will be addressed in the future study. 
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FIG. 1. (a) We consider a discrete folding of the triangular lattice. The fold angle (with respect 
to the adjacent triangular plaquettes) is discretized into either 6 = or n. (b) A drawing of a 
transfer- matrix strip is shown. 
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FIG. 2. The free-energy gap ([7]) is plotted for the bending rigidity K and the system sizes 6 < L < 
14. 
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FIG. 3. The transition point K(L) ([8]) is plotted for 1/L 2 . The linear least-squares fit for (+) 
L = and (x) 1,2 mod 3 (6 < L < 14) yields K = -0.2697(12) and -0.2695(14), respectively. 
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FIG. 4. The latent heat Q(L) §U§ is plotted for 1/L 2 . The linear least-squares fit for (+) L = 
and (x) 1,2 mod 3 (6 < L < 14) yields Q = 0.0482(59) and 0.0391(38), respectively. 
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FIG. 5. The free-energy gap ([7]) is plotted for the bending rigidity K and the system sizes 6 < L < 
14. Tentatively, the defect parameter ([5]) is set to K' = 0. 
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